Magnetic order, magnetic correlations and spin dynamics in the pyrochlore 

antiferromagnet Er 2 Ti 2 Oy 
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Er2TiaOr is believed to be a realization of an XY antiferromagnet on a frustrated lattice of 
corner-sharing regular tetrahedra. It is presented as an example of the order-by-disorder mech- 
anism in which fluctuations lift the degeneracy of the ground state, leading to an ordered state. 
Here we report detailed measurements of the low temperature magnetic properties of Er2Ti207, 
which displays a second-order phase transition at Tn — 1.2 K with coexisting short- and long-range 
orders. Magnetic-susceptibility studies show that there is no spin-glass-like irreversible effect. Heat- 
capacity measurements reveal that the paramagnetic critical exponent is typical of a 3-dimensional 
XY magnet while the low-temperature specific heat sets an upper limit on the possible spin-gap 
value and provides an estimate for the spin-wave velocity. Muon spin relaxation measurements 
show the presence of spin dynamics in the nanosecond time scale down to 21 mK. This time range 
is intermediate between the shorter time characterizing the spin dynamics in Tb2Sn207, which also 
displays long- and short-range magnetic order, and the time scale typical of conventional magnets. 
Hence the ground state is characterized by exotic spin dynamics. We determine the parameters of a 
symmetry-dictated Hamiltonian restricted to the spins in a tetrahedron, by fitting the paramagnetic 
diffuse neutron scattering intensity for two reciprocal lattice planes. These data are recorded in a 
temperature region where the assumption that the correlations are limited to nearest neighbors is 
fair. 

PACS numbers: 75.40.-s, 75.10.Dg, 76.75.+i 



I. INTRODUCTION 



Because of the geometrical frustration of their mag- 
netic superexchange interactions, the insulating py- 
rochlore compounds R2M2O7, where R stands for a 
magnetic rare earth ion and M = Ti or Sn, display a 
variety of unusual magnetic behaviorsj^ Examples in- 
clude (i) spin ice systems Ho2Ti 2 07 and Dy 2 Ti 2 07^ 
(ii) Yb 2 Ti 2 07 with a sharp transition in the spin dy- 
namics finger-printed by a pronounced peak in the spe- 
cific heat; 4 * and (hi) Tb 2 Sn 2 07 in which magnetic Bragg 
reflections are observed at low temperature by neutron 
diffraction^ 5 - while no spontaneous magnetic field is found 
by the muon spin rotation (/LtSR) technique.- In addition, 
even when a spontaneous field and magnetic Bragg re- 
flections are detected, as it is expected for a conventional 
ordered magnet, persistent spin dynamics in the ordered 
state are surprisingly observed, e.g. in Gd 2 Ti 2 07 and 
Gd 2 Sn 2 07. 7 ~ 13 In terms of crystal-field anisotropy, the 
spin-ice systems are strongly Ising-like. Tb 2 Sn 2 07 has 
also an Ising anisotropy, but not so strong. Yb 2 Ti 2 07 is 
XY-like from the crystal-field point of view and the Gd 
compounds are approximately isotropic. 

Although Yb 2 Ti 2 07 is XY-like, its magnetic moments 
are not perpendicular to the local <111> axes and it does 
not display long-range magnetic order^i 4 - although this 
absence of order has been dispute d 15 ! 16 and is still un- 



der debated These anisotropy and absence of long-range 
order also pertain for Yb 2 GaSb07»i& Hence, it was of 
great interest when Er 2 Ti 2 07 was reported to be XY-like 
and to display long-range order at low temperature with 
the Er 3+ magnetic moments perpendicular to their lo- 
cal [111] axes^ Later on, however, coexisting short- and 
long-range orders were found and soft collective modes 
were detected^ The presence of the soft modes has been 
attributed to the incommensurate value of the canting 
angle. 21 These astonishing inferences call for more de- 
tailed data and analysis. This is the purpose of this work. 
One of the experimental advantages of Er 2 Ti 2 07 over 
Yb 2 Ti 2 07 and Yb 2 GaSb07, is the possibility to produce 
large high-quality crystals. 

Another reason for the interest in the Er 2 Ti 2 07 system 
is the following. As a realization of an XY antiferromag- 
net on a pyrochlore lattice, it is a natural candidate for 
observing the phenomenon of order by disorder which has 
been discussed theoretically for more than three decades 
since the pioneering study by Villain and coworkers. 22 
Bramwell et al. indeed showed that while the zero tem- 
perature ground state is degenerate, thermal fluctuations 
select a subset of the manifold and induce a first order 
phase transition to a conventional Neel ground stated 
The order by disorder mechanism has been confirmed in 
several subsequent works, see e.g. Refs. [l9|, [2J-|26|, and 
interestingly the more recent studies consider the effect 
of quantum fluctuations and tend to explain the second 
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FIG. 1. (color online). The network of corner sharing regular 
tetrahedra formed by the rare earth atoms in the pyrochlore 
structure in which Er2Ti207 crystallizes. The axis of trigonal 
symmetry at the position of a rare earth is one of the cube 
diagonals. There are two types of tetrahedra in the network, 
which differ by their orientation: type B is rotated by 90° 
about the cubic axes with respect to type A. We distinguish 
the two sets by two colors in the drawing. Since each rare 
earth is at a corner shared by two tetrahedra, one of each 
kind, either the set of the four corners of all the A tetrahedra 
or the set of the four corners of all the B tetrahedra is sufficient 
to describe the Er 3+ lattice. 



order nature of the transition experimentally observed in 
Er 2 Ti 2 7 . 

The organization of this paper is as follows. Section [TT1 
gives a survey of the physical properties of Ei - 2Ti207 and 
discusses its magnetic structure. In Sec. UTT] we describe 
the growth of the single crystals, their basic characteriza- 
tions and the experimental methods used in the present 
work. Section IIVI presents our investigation of the bulk 
properties of the system, including magnetic suscepti- 
bility and specific heat measurements and their analysis. 
The following section (Sec. [V} deals with the microscopic 
techniques, i.e. muon spin relaxation and neutron scat- 
tering in the paramagnetic phase. A summary of our 
key results is given in Sec. I VII The physics of effective 
one-half spins (Kramers doublets) on a tetrahedron is 
described in Appendix [5] Appendix [B] outlines the cal- 
culation needed for the analysis of the specific heat data 
presented in Sec. |lVl 



II. PHYSICAL PROPERTIES OF Er 2 Ti 2 7 AND 
MAGNETIC STRUCTURE 

Erbium titanate, Er2Ti207, is an insulating pyrochlore 
compound that crystallizes into the cubic space group 
Fd3m, with the lattice parameter a = 10.0727 (1) A at 
room temperature and x = 0.3278(8), the free position 
parameter allowed by the space group which character- 
izes the 48f site occupied by oxygen^ 

The Er 3+ ions which occupy the 16d Wyckoff posi- 
tions in the space group, are located at the vertices of a 
corner-sharing network of tetrahedra; see Fig. [T] A sin- 
gle tetrahedron with four Er sites comprises the primitive 
unit cell. The Er 3+ crystal sites are all equivalent and 
the local symmetry is Z?3d, where the 3- fold axes pass 



through the center of a tetrahedron, in the directions 
[111], [111], [111] and [III] for the four sites numbered 
1, 2, 3, and 4, respectively on a single tetrahedron. Ac- 
cording to Hund's rules, the total angular momentum of 
the Er 3+ ion in its ground multiplet is J = 15/2. The 
16-fold degeneracy is lifted into Kramers doublets by the 
crystal electric field (CEF). 

The ground state doublet can be described as an ef- 
fective spin 5 = 1/2 with a quantization axis z paral- 
lel to the local trigonal axis. It is well isolated from 
the excited doublets, the lowest being at about 74 K 
above the ground state in temperature units We shall 
write the CEF ground state doublet wavefunctions as 
\4>^). This doublet is characterized by its spectroscopic 
factors along and perpendicular to the trigonal axis, 
gu and g± respectively. From a global analysis of the 
CEF for the pyrochlore -R2T12O7 series it has been de- 
duced that gu = 2gj\((pQ\Jz\<f*o}\ — 1-8(5) and g± = 
gj\(<Po \ J±\<t>a )\ = 7.7 (1)M Here gj = 6/5 is the Lande 
factor. These spectroscopic factors are related to matrix 
elements that we shall need for the analysis of neutron 
scattering data. We have 

jcef = = -<^»o Wo) = 0.75(20), (1) 

t CEF = (4\J±\cft) =6.42(8). (2) 

where J± = J x ±iJ y . The matrix elements obviously refer 
to quantities written in local axes; see Appendix IA II for 
a discussion. By definition, the large difference between 
Jcef and icEF (and obviously also between g^ and g±) 
reflects the strong CEF anisotropy of the XY type of the 
Er spins (in contrast to the Ising limit for which j'cef 3> 
^CEf)- 

The compound displays a magnetic phase transition 
at Tn — 1.2 K32 The large negative value of the Curie- 
Weiss temperature 8cw (#cw = —22 K is deduced from 
susceptibility data measured between 20 and 50 K; see 
Refs.H^ and HO) suggests a strong antiferromagnet cou- 
pling. 

Neutron diffraction shows the magnetic structure be- 
low Tn to be noncollinear with the propagation vec- 
tor k = (0, 0,0)ji^ From polarized neutron diffraction, 
the Er 3+ magnetic moment is determined to be m = 
3.25 (9) /iB at low temperature^ A note of caution seems 
justified at this juncture: m is not directly related to 
the spectroscopic factors which have been determined for 
a paramagnetic ion, since the molecular field has to be 
taken into account for an estimation of m. The diffrac- 
tion data have been originally describe d 19 i 31 with the 
irreducible representation^ We notice that this descrip- 
tion also considered in Refs. [2^, [3^, andHH has been re- 
cently disputed by Briffa et al£^ 

In fact the available microscopic information provides 
an insight into the moment orientation. Let us consider a 
one-half spin subjected to a molecular field oriented at a 
polar angle 9 from the [111] axis. The following relation 



3 



100 



£ 10 



o 
O 



Er 2 Ti 2 7 

50 mK, zero field 
(2,2,0) peak 




Instrument 
resolution 




0.001 



0.01 

k - 9(2,2,0)1 (A" 



0.1 



CO 



sz 

CD 



80000 
60000 
40000 
20000 




1 1 1 1 1 1 

— — background 

— calculated intensity 


ii' 

Er 2 Ti 2 7 


• observed intensity 


— difference 




- 




(222 








. — 






Cm I? O rj i 
T VO hf O (N _ 




o d 


„S^i|S SS 22-Jo 


: 1 \ 




GIG 1 pSp ' ps - 




lit 


i 1. iirr.r iir a 




_H ^ ~ 

1,1,1,1,1,1,1 



20 40 60 80 100 
Scattering angle 29 (deg.) 



120 



FIG. 2. (color online) Neutron intensity for the (2,2,0) Bragg 
peak of Er2Ti207 plotted as a function of \q — 5(2,2, o)|- The 
data are from Ruff et aZ.— The instrumental resolution is de- 
termined from the reflection measured under a 3 T magnetic 
field. The line for \q — 17(2,2.0) I > 0.015 A -1 is described in the 
main text. 
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tan 9 



4m 2 



4rn 2 - [i\g\ ' 



(3) 



Numerically this gives 8 « 20° . This means that the field 
is not far from being parallel to the [111] axis. However, 
the polar angle of m is given by tan <j) = (g±/g\\ ) 2 tan 9, 
i.e. <f) ~ 80°. Taking into account the uncertainties on 
g±, (/ii, and m, this analysis indicates that the moment 
is perpendicular to the local [111] axis, or at least close 
to being perpendicular. This is consistent with the mag- 
netic structure first proposed b y C hampion et al.^ and 
invalidates the proposal of Ref. l2ll . 

Regardless of the magnetic structure, the (2,2,0) Bragg 
reflection has been shown to be anomalous with broad 
tails. 20 Quantitatively, we find that the intensity around 
this position is proportional to I?— <j(2,2,o) I with rj ~ 0.6 
for \q - ^(2,2,0)1 > °- 03 A -1 ; see Fig. [2] The value of the 
exponent rj is much reduced compared to rj ~ 1.8 found 
in Tb2Sn207.— ^ For a conventional ordered compound 
the neutron intensity would be Gaussian-like, i.e. with a 
much steeper slope. 

Inelastic neutron scattering data recorded in zero field 
suggest the presence of magnetic soft modest in agree- 
ment with the power law behavior of the specific heat 
below Tn-~ At first sight this is surprising given the ex- 
pected strong crystal-field anisotropy of the Er 3+ ions. 21 



III. EXPERIMENTAL 

Er2Ti2C>7 single crystals were grown by the float- 
ing zone technique using a commercial optical furnace. 
Feed rods were prepared from high purity oxides (Ti02, 
99.995% and Er 2 3 , 99.99%), mixed and heat treated up 
to 1180°C. After sintering, a rod was heat treated up to 



FIG. 3. (color online). An example of an x-ray diffraction 
pattern recorded at room temperature for an Er2Ti207 pow- 
der obtained after crushing part of a crystal. The radiation 
used is Co K Q . The line at the bottom of the graph shows 
the difference between the data and the refinement model. 
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FIG. 4. (color online). Specific heat versus temperature for 
three samples cut from an Er2Ti2 07 crystal prepared as in- 
dicated in the main text. For the first sample, no subsequent 
treatment has been performed, while the other two result from 
post-growth heat treatments as follows. Case #1: heat treat- 
ment for 7 days under oxygen at 1150°C, followed by a slow 
cooling down to 400° C. Case #2: heat treatment for 7 days 
under argon at 1150°C, followed by a slow cooling down to 
400°C. 



1350°C in air. Crystal growth conditions were optimized 
under air (1 !/mn) at the growth rate of 2 mm/h coupled 
with a rotation rate of 30 rounds per minute. As with 
most of the titanate pyrochlores, Er 2 Ti 2 07 crystallized 
rods are mostly transparent, with a slight pink color. No 
phases other than the cubic one with Fd3m space group 
were detected by x-ray powder diffraction experiments; 
see Fig. [3] for an example. As grown and post growth 
heat-treated crystals were characterized by specific heat 
measurements. No noticeable differences were detected; 
see Fig.Hl This is in contrast to the case of Tb 2 Ti 2 07j2£ 
The position of the peak provides a measure of the critical 
temperature. We get Tn = 1-23(1) K. For comparison 
the following values have already been published: 1.25 K 
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(Ref. [H) and 1.173 (2) K (Ref. QJ) in reasonable agree- 
ment. 

The investigation of the macroscopic properties of the 
system consisted of magnetic susceptibility experiments 
performed with a commercial magnetometer (Magnetic 
Property Measurement System, Quantum Design inc.) 
down to 2 K, and of heat capacity measurements. For this 
latter physical property, the temperature range from 0.48 
to 20 K was investigated with a commercial calorime- 
ter (Physical Property Measurement System, (PPMS) 
Quantum Design inc.) equipped with a 3 He stage us- 
ing a standard thermal relaxation method. Additional 
measurements between 0.11 and 2.50 K were performed 
with a home-made calorimeter inserted in a dilution re- 
frigerator using a semi-adiabatic technique. 

The fiSK measurements were carried out at the Euro- 
pean Muon Spectrometer of the ISIS facility (Rutherford 
Appleton Laboratory, United Kingdom) and the Low 
Temperature Facility of the Swiss Muon Source (S/xS, 
Paul Scherrer Institute, Switzerland). The muon beam 
is pulsed at the former facility and pseudo-continuous at 
the latter. 

The neutron scattering experiments were performed 
at the Institut Laue Langevin (ILL, Grenoble) with the 
lifting-counter diffractometer D23 of the CEA collabo- 
rating Research Group (CRG). 



IV. BULK MEASUREMENTS 

Here we shall first discuss magnetic susceptibility mea- 
surements. Then we shall present zero-field specific heat 
results and finish with the determination of the magnetic 
phase diagram using specific heat-data recorded under 
magnetic fields. 



A. Magnetic susceptibility 

Since the magnetization measurements were performed 
on a needle-shape sample and the field was applied along 
its long axis, the demagnetization field is negligible. 

Classically, the static magnetic susceptibility \ is ex- 
pected to follow a Curie- Weiss law far from the ordering 
temperature in the paramagnetic regime. It reads 



C 



X. 



T — 8qw 



(4) 



where the Curie constant C can be expressed in terms of 
the so-called paramagnetic moment m para : 
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I Horn 



para 



3k T 



(5) 



where v = a 3 /N cc \i with N ce \\ being the number of Er 3+ 
ions in the cubic cell (iV ce ii = 16). For an isolated Er 3+ 
ion, m para = g.j^J{J + 1)mb = 9.58 jLt B - 
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FIG. 5. (color online) Inverse of the magnetic susceptibility 
of an Er2Ti207 crystal versus temperature in a large temper- 
ature range. The solid line results from a fit with the Curie- 
Weiss law. A field of l mT is applied along a diagonal of the 
cubic crystal structure. The inset shows the low temperature 
range of the data. 



In Fig. [5] we display our result for the inverse of the 
static susceptibility versus temperature in a large tem- 
perature range. The Curie- Weiss law provides a good 
description of our data above 30 K. The fit gives for 
the Curie- Weiss temperature #cw = —17.5(3) K and 
C = 3.73 (4) K. This means that m para = 9.55 (10) ^b, 
in agreement with the result for an isolated Er 3+ ion. 
Because 9qw is negative, the dominant exchange inter- 
actions are antifcrromagnetic. For comparison, fitting 
data recorded between 20 and 50 K on a powder sam- 
ple in a field of 1 mT, Bramwell et at 30 found val- 
ues of 6>cw = —22.3(3) K and m para = 9.34 (9) /is, 
For data recorded between 50 and 300 K in a field of 
50 mT these authors obtain 6cw = —15.93(3) K and 
Wpara = 8.936 (4) /1b, the latter at variance with the re- 
sult for an isolated Er 3+ ion. From our result, we com- 
pute for the frustration index^ 1 / = |#cw|/7n = 14. 
Since / ^> 1, Ei'2Ti207 is a strongly frustrated mag- 
net. In addition, assuming the Er 3+ magnetic moments 
to interact through a simple nearest-neighbor Heisenberg 
interaction with exchange integral I (I > 0), i.e. 



n = \ £ J > 



the molecular-field approximation predicts 
„ 3/c B |6»cw| 



z Dn J{J+ 1)' 



(6) 



(7) 



We denote as z nn the number of nearest neighbor Er 3+ 
ions to a given Er 3+ ion. In our case z nn = 6. From 
the measured 6cw value and taking into account that 
,/ = 15/2, we compute l/k B = 0.138 (2) K. 

We have also measured the susceptibility for 2.0 < T < 
6.0 K under a field of 1 mT applied along a [111] axis 
using two protocols; see Fig. [51 Contrary to a previ- 
ous report^ we do not observe any history dependent 
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FIG. 6. (color online) The magnetic susceptibility of an 
Er2Ti2C>7 crystal versus temperature in the low temperature 
region. The same magnetic response is observed using either 
the field or zero-field cooling procedure, fc and zfc respec- 
tively. A field of f mT is applied along a diagonal of the cubic 
crystal structure. In the zfc mode the sample was warmed to 
300 K and then cooled to 2 K before applying the magnetic 
field. In the fc mode the field was applied at 300 K and the 
sample was subsequently cooled down. 

effect at T < 3.2 K. Hence, there is no spin-glass-like 
irreversible effect for our Er 2 Ti 2 07 crystals. 
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FIG. 7. (color online) Low temperature specific heat of 
Er2Ti2 07. The open symbols show our experimental data 
compared to literature results from Blote et al^ and Sid- 
dharthan et aZ4£ Note the large overlap between our results 
obtained with the PPMS and the home-made calorimeter in- 
serted in a dilution refrigerator. The filled symbols present the 
electronic specific heat deduced from the data of Siddharthan 
et al. and this work. The dashed lines show the contributions 
from the nuclear specific heat and the full line results from 
a fit of the low temperature electronic specific heat to the 
Ceioc = BT 3 law, with B = 2.50(15) JK _4 mor 1 . The in- 
set displays the very low temperature details. The solid lines 
result from fits as explained in the main text. 



B. Specific heat in zero magnetic field 

Here we present and discuss zero-field specific heat 
data recorded for Er 2 Ti 2 C>7 crystals. It is well known 
that they may lead to a characterization of the low energy 
magnetic modes, detect indirectly a dynamical magnetic 
component in the ordered state, determine the universal- 
ity class of the system under study, and gauge a possible 
residual entropy at low temperature. 

In Fig. [7] we display our data, in the low temperature 
regime. While our results are in reasonable agreement 
with data published by Siddharthan et al., 40 they differ 
drastically at low temperature with the ones of Blote 
et al£2- We do not understand the origin of the large 
difference with the data of Blote et al. In the following 
we shall focus on the analysis of our data and the ones of 
Siddharthan et al. Note that here we have not considered 
the results of Sosin et al. 34 since they do not extend to 
very low temperature. 

In the low temperature range shown in Fig. [7] the spe- 
cific heat arises from three origins: the contribution from 
the nuclei with non-zero spins at low temperature, then 
the low energy magnon modes at higher temperature and 
the critical fluctuations around Tn- As a start we shall 
focus on the first two contributions, i.e. the specific heat 
of nuclear and magnon origins, C n and C sw , respectively. 

We begin with providing a theoretical background 
for these contributions, focusing first on the nuclear 
one. C n arises from the 167 Er nuclear magnetic mo- 
ments ( 167 Er is the only non-spinless isotope of Er 



with 23% natural abundance), since the contribution 
of the two Ti isotopes is negligible, as in the case of 
Tb 2 Ti 2 7 . 38 Contrary to Tb 2 Ti 2 7 the quadrupole in- 
teraction is not negligible compared to the Zeeman in- 
teraction. This is due to the fact that the quadrupole 
moment Q W7 of 167 Er is larger than that of 159 Tb (3.565 
vs 1.432 barns) and the gyromagnetic ratio jier of 167 Er 
is much smaller, in absolute value, than that of 159 Tb 
(-7.7157 vs 64.31 Mrads" 1 ! 1 - 1 ); see Ref . [H The Zee- 
man and quadrupolar Hamiltonians are written 

■Hzcc = — S-7167I • Bhyp (8) 

and 

H Q = hu Q [3I 2 z -1(1+1)], (9) 

respectively. In these equations, I is the 167 Er spin op- 
erator (I = 7/2) and Hujq = if^i-T) wnere V xx is the 
principal component of the electric field gradient tensor 
acting on the rare earth nucleus with z being as before 
the local three-fold axis. The symmetry at the rare earth 
site imposes the electric-field gradient to be axial. Be- 
cause the Er 3+ ordered magnetic moments are (nearly) 
perpendicular to z we shall also take Bh yp perpendicu- 
lar to this axis. As usual fi and e stand for the Dirac 
constant and the proton electric charge, respectively. 

The nuclear energy levels are determined after diago- 
nalization of the Hamiltonian H n = "Hzcc + ^q and C n is 
readily obtained. "H n depends on two parameters, B^ yp 
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and V zz . While an estimate for V zz is provided in Ap- 
pendix [B] -Bhyp w ill be a fitting parameter. 

The other contribution to the low temperature specific 
heat arises from magnons. Low energy magnons have in- 
deed been observed in neutron scattering experiments. 20 
The dispersion relation 7kj(q) for their lowest energy 
branch is needed to compute C sw . An approximate ex- 
pression valid at small wavevectors is 

&V(q) = &V(?) = A 2 SW + h 2 v 2 sw q 2 . (10) 

Here A sw is the gap energy of the magnon spectrum at 
the zone center and v svr is the magnon velocity. We 
note that a dispersion relation has recently been pro- 
posed for Er2Ti2C>7 in the framework of linear spin- wave 
theory— The applicability of this theory in frustrated 
systems might be questionable as recently discussed in 
the case of the triangular lattice. 42 Still, the model of 
Ref. leads to an anisotropic dispersion relation. The 
resulting specific heat depends on a single magnon veloc- 
ity which is the geometrical mean of the three magnon 
velocities along orthogonal axes. In our model it corre- 
sponds to v sw . 

When A sw is negligible, the magnon specific heat C sw 
can be computed in the temperature range where only 
small wavevectors are at play, i.e. when Eq. [10] applies. 
The expected T 3 law for C sw is derived: 

C sw = AT 3 with A = ^A/- A -3^, (11) 

where A/a is Avogadro's constant. This result only holds 
at sufficiently low temperature. 

Having established the theoretical background, we now 
perform the specific heat data analysis. We shall do it in 
two steps. 

We first attempt to determine whether a T 3 behavior 
can be observed. A fit to the measured specific heat at 
the lowest temperatures for which the magnon contribu- 
tion should be negligible enables us to estimate C n (T) 
and then to subtract it from the measured heat capacity. 
The resulting electronic heat capacity, C c i C c, is presented 
in the main panel of Fig. [7J for our data and the ones 
of Siddharthan et al. It follows nicely a T 3 law in a re- 
stricted temperature range, but deviates above w Tn/2.5, 
in contrast to published resultsj 19 i 34 This observation jus- 
tifies to identify C c i cc with C sw . The T 3 behavior is not 
expected to be seen at low temperature if the energy gap 
is appreciable. The effect of the gap might be seen around 
T = 0.2 K; see Fig. [JJ However, C c i ec becomes very small 
at that temperature and difficult to measure as reflected 
by the distribution of the C c i cc data. Hence, we cannot 
determine whether a gap is present from this plot. Nu- 
merically, since we can identify B given in the caption of 
Fig. [7J with A of Eq.ELIl we get v sw = 86(2) ms" 1 . 

The second step for the interpretation of the specific 
heat data consists in fitting the measured specific heat to 
the sum C n + C sw . This sum depends on three param- 
eters, A sw , v sw and £?hyp- The fit is shown in the inset 
of Fig. [JJ Its temperature range is restricted on the high 



temperature side because of the requirement that only 
the low energy magnons determine the value of the inte- 
gral. Two solid lines are drawn in the figure since the two 
data sets are slightly different at low temperature. Both 
sets can be fit to a range of gaps extending from nearly 
zero to an upper bound. We find A sw //cb < 0.5 (1) K for 
both data sets. This is consistent with the value of Sosin 
et al£^ We also derive v sw = 84 (2) and 82 (2) ms" 1 and 
B hyp = 345 (10) and 305 (5) T for the Siddharthan et al. 
data and our data, respectively. Note that -Bhyp depends 
very little on the actual value chosen for V zz . 

We now discuss these results, starting with the bound 
on the spin gap energy. This bound is really small and 
might be surprising at a first sight given the strong mag- 
netic anisotropy of Er2Ti2C>7. However, the Er 3+ mag- 
netic moment lies at a polar angle of nearly 90° with re- 
spect to the local three-fold axis. This angle is imposed 
by the relatively strong crystal field interaction. There 
is still a continuous degree of freedom for the azimuthal 
angle. The magnetic order breaks this rotational sym- 
metry and A sw is a measure of the residual anisotropy 
energy. We note that the upper bound value for A sw is in 
the expected range if it arises from the dipole interaction 
between the Er 3+ magnetic moments. 

We examine now the magnon velocity value and tenta- 
tively relate it with the exchange integral introduced in 
Eq. |51 For this purpose we resort to the phenomenologi- 
cal dispersion relation 

ftV(q) = h 2 u 2 {q) = A S 2 W + [Xz an Jsm{qd)] 2 , (12) 

where d is the distance between two magnetic atoms. In 
fact, had we written sm(q z d) instead of sm(qd), Eq. Q2] 
would give the dispersion relation of an antiferromag- 
netic chain running along the Z direction and of lattice 
parameter d for which the number of nearest neighbors 
is z nn = 2; see for example Ref. H^. Here we shall take d 
as the distance between two Er atoms, i.e. d = a/ (2v2). 
Identifying the small q expansion of Eq. [T2] with Eq. [TU] 
we have 

<W = ^a~, (13) 

using J = 15/2. This relation leads to I/fcs = 
0.118 (4) K, a value in reasonable agreement with the one 
derived from the Curie- Weiss constant; see Sec. II V Al We 
shall return to the interpretation of the magnon velocity 
at the end of Sec. EH once a nearest neighbor Hamilto- 
nian consistent with the symmetry of the lattice has been 
introduced. 

It is possible to get interesting information from the 
-Bhyp values. Using the hyperfine constant which is 
87 (1) T/z-tBr 4 - the magnetic moment at the origin at the 
field can be obtained. We derive for the moment 4.0 (2) 
and 3.5 (1) /ie from the Siddharthan et al. measurement 
and ours, respectively. The latter value is consistent 
with the neutron result^ Hence, contrary to Tl^S^C^ 
(Ref. |45| ) there is no reduction of the hyperfine field at the 
167 Er nuclei due to electronic spin dynamics. This means 



7 




FIG. 8. (color online) Zero-field specific heat of Er2Ti2C>7 
versus the reduced temperature parameter r in the param- 
agnetic regime for three data sets. The full (dotted) line is 
the prediction of Eg. [Til for a = -0.015 (-0.134). For a = 
-0.015 we find C sh = 1.7 (1) JK" 1 moP 1 . The critical regime 
is observed up to r ~ 0.2. 



that the characteristic time for the electronic spin-flip is 
substantially larger than the 167 Er spin-lattice relaxation 
timeAi 

We now turn our attention to the critical behav- 
ior of the specific heat in the paramagnetic phase. In 
Fig. [5] we display our data using a reduced temperature 
scale. We expect to observe the usual power law critical 
behavior: 46 ' 47 



C e \ec(T) 



Csh 



- 1 



(14) 



where C s h is a constant and a the specific heat critical 
exponent. By definition, C e \ ec (T) has a maximum at Tn. 
This enables us to determine Tn, as already mentioned 
in Sec. IIIIl The exponent a is expected to be a = —0.015 
and —0.134 for the three-dimensional XY and Heisenberg 
magnets, respectively. As seen in Fig. [5J Eq . IT41 provides 
a better account of the data for the XY expected. 

Before leaving this section, we discuss the entropy 
variation of our system using our specific heat results. 
We have extended the specific heat measurement C p of 
Er2Ti2C>7 up to approximately 20 K and measured the 
specific heat of the isostructural non-magnetic compound 
Y2Ti2C>7 in the same temperature range. The sum of 
the contributions from the nuclear moments and the lat- 
tice, the latter being estimated from scaling the Y2Ti2C>7 
result,— is presented in the inset of Fig. [HI as well as 
C P (T). The resulting C i ec (T) obtained from subtracting 
the contributions of the nuclear moments and the lattice 
to C P (T) is displayed in the main frame of Fig. [5] In 
Fig. [JTJwe present AS c \ ec (T), which is the temperature 
variation of the electronic entropy obtained by integrat- 
ing C e iec(T')/T' from 0.115 K, the lowest measured tem- 
perature, to the temperature of interest T. Remarkably, 
AS e \ ec (T) reaches the i?ln(2) value at approximately 8 K. 
This is the entropy variation expected for an isolated dou- 



FIG. 9. (color online) The electronic specific heat of Er2Ti2C>7 
in an extended temperature range. In the inset are displayed 
the total specific heat (C p ) of Er2Ti2 07 and the sum of the 
estimated nuclear and lattice contributions. 
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FIG. 10. (color online) The variation of the electronic entropy 
AS'eiec(T) of Er2Ti2 07 versus temperature T. The solid line is 
computed assuming two singlet states at energies correspond- 
ing to and 2.85 K and two doublets at 74 and 85.8 K. A 
signature of the two doublets has been seen by inelastic neu- 
tron scattering.— The ground state doublet has been split to 
account for the magnetic order of Er2Ti2 07 below Tn, with 
the splitting taken as a fitting parameter. Admittedly, this 
is a very rough description of the physics. The dashed line 
is the result of the computation of AS e \ ec (T) when only the 
first two levels are taken into account. 



blet state. Therefore the residual entropy left as T — » 
is extremely small if any. For temperatures above 8 K, 
AS , e i ec (T) keeps increasing owing to the contribution of 
the excited CEF levels. This is clearly shown in Fig. [TUl 
where a comparison is made between the results of com- 
putations of AS'cioc (T) either taking into account the first 
excited CEF levels or neglecting them. 



C. Specific heat under an external magnetic field 

We have constructed the phase diagram of Er 2 Ti 2 07 
in the field-temperature plane using specific heat data. 
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FIG. 11. (color online) Temperature dependence of the spe- 
cific heat of a Er2Ti2C>7 single crystal for different magnetic 
field intensities applied along [110]. The maximum of the 
specific heat peak moves to lower temperatures as the field 
increases up to 1.7 T. No peak is observed when the field 
strength is above 1.7 T. 
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FIG. 12. (color online) The phase diagram derived from spe- 
cific heat measurements for the three main crystal directions 
of cubic Er2Ti207. The dashed-dotted lines are guides to the 
eye. 



Examples of measurements are shown in Fig. 1111 For a 
given external field we have determined the temperature 
at which the specific heat displays a maximum. The po- 
sition of the maximum as a function of the field intensity 
for a given field orientation relative to the crystal axes is 
displayed in Fig. Q21 Our results are consistent with the 
ones already published , 20 i 34 but only qualitatively. Note 
that here we establish the phase diagram for the three 
main directions of a cubic compound. These data sug- 
gest a quantum critical point to be present slightly above 
2 T and to be dependent on the field orientation. 



MICROSCOPIC TECHNIQUE 
MEASUREMENTS 



Single crystals of Er2Ti2C>7 have been studied by two 
microscopic experimental techniques: positive muon spin 
relaxation (//SR) and neutron scattering. We shall first 
discuss the /j,SR results. 



A. /iSR 

In the longitudinal geometry that we have used, a /iSR 
spectrum recorded in the magnetically ordered state of 
a crystal is expected to display either (i) at least one 
damped oscillation if the initial muon beam polariza- 
tion is not parallel to the spontaneous field at the muon 
site, or (ii) a missing fraction if the oscillation cannot be 
resolved42 None of these two possibilities was observed at 
ISIS or S/iS. The zero-field spectrum recorded at 21 mK 
is displayed in Fig. the spectral shape shows little 
change up to w 0.5 K. The same type of spectra was 
also observed for different orientations of the initial muon 
beam polarization relative to the crystal axes. Hence, the 
absence of oscillation cannot be attributed to the initial 
muon beam polarization which would be parallel to the 
internal field. This situation would moreover be unex- 
pected because of the four equivalent symmetry (111) 
axes at the Er 3+ site. In addition to the absence of oscil- 
lation, the shape of the zero- field spectra is extremely un- 
usual. Indeed, the spectral slope is quasi-constant up to 
~ 0.4 /is, then it increases and eventually monotonically 
decreases above « 0.7 /is. This kind of behavior drasti- 
cally differs from the usual relaxation spectra, which are 
characterized by a monotonically decreasing slope, as e.g. 
in the common exponential relaxation. It does not either 
remind the shape of spectra associated with a static or 
quasistatic field distribution at the muon. At this stage 
we cannot actually conclude whether the spectral shape 
is mainly influenced by a static field distribution or dy- 
namical effects. We shall come back to this point when 
commenting the data recorded in applied fields. 

Since no spontaneous muon precession is observed, it is 
tempting to make an analogy with Tb2Sn2 07 for which 
there is also no detected oscillation in the magnetically 
ordered stated This analogy would be completely justi- 
fied if the zero-field spectrum was exponential. However, 
we have just remarked that it is obviously not the case. It 
is not Gaussian either, as first suggested^ A close look 
at the data in Ref. shows that the zero-field relaxation 
is in fact consistent with the one we observe. 

Before discussing further the significance of our re- 
sult, it is worthwhile to consider the longitudinal field 
spectra shown in Fig. Q21 First note that with the con- 
sidered field intensities the compound has a priori not 
crossed the phase diagram boundary displayed in Fig. [121 
While the application of fields below 0.5 T has little in- 
fluence on the spectral shape, it is not the case for higher 
fields. Interestingly, for the highest fields shown in the 
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FIG. 13. (color online) ^iSR spectra of a Er2Ti2 07 crystal 
taken at S/iS with the initial muon beam polarization parallel 
to the [111] crystal direction. The spectra were recorded in 
zero and various longitudinal fields as indicated in the figure. 



in Er2Ti2 07 are slower by roughly an order of magnitude. 
From neutron spin-echo measurements it is known that 
in the ordered state of Tb2Sn2 07 spin correlations near 
the zone center are static^ within the technique time 
scale, while they are dynamical in nature far outside the 
center of the zone^ It would be worthwhile to examine 
the dynamics of the magnetic correlations in Er2Ti207 
with the neutron spin-echo technique. 

From the analysis of the nuclear specific heat in 
Sec. IIVB1 it was inferred that the ratio of Tf to the nu- 
clear spin-lattice relaxation time was larger in Er2Ti207 
than in Tb2Sn207. Assuming the spin-lattice relaxation 
times in the two compounds to be similar, the nuclear 
specific heat data are consistent with the results of the 
analysis of the /iSR spectra. 



B. Neutron scattering 



figure, the spectra tend to be described by an exponen- 
tial function; 51 This implies that the muon repolarization 
is not at the origin of the field dependence of the spec- 
tra and that dynamics is mainly influencing the muon 
response in Er2Ti207. The time scale of this dynam- 
ics can be roughly estimated. For this purpose we have 
recourse to the Lorentzian field dependence of the expo- 
nential relaxation rate Xz in the motional narrowing, i.e. 
fast fluctuation limit i 49 > 52 In this model the relaxation 
rate in a field B cxt is such that Xz(B cxt )/ Xz(B cxt = 0) 
— 1/2 for B cxt — B X j 2 — 1/(7^ Tf). Here Tf is the fluctu- 
ation time of the spins and 7^ = 851.615 Mrads -1 T -1 
is the muon gyromagnetic ratio. Taking 1 T as an order 
of magnitude for -B1/2 we find Tf w 10~ 9 s. 

We have therefore established that the muons are prob- 
ing dynamical fields. No oscillation is detected in zero- 
field because the mean field at the muon site does not 
keep a constant value for a time sufficiently long for a 
muon spin precession to be observed^ We denote this 
field as Bfl uc - It arises from the dipole interaction of 
the muon magnetic moment with the Er 3+ magnetic mo- 
ments. Given the size of the Er 3+ magnetic moment, we 
estimate -Bfl uc in the range 0.1-0.2 T. 

With our knowledge for Tf and -Bfluc, we find 
T i 7m Bfl uc = -Bfluc/-B cx t <C 1. For the simple model of 
Bfluc flipping from parallel to antiparallel to an axis per- 
pendicular to initial muon polarization S^, we would 
expect the zero-field relaxation to be exponential^ 9 - as 
found for Tb2Sn2C>7. However, experimentally this is not 
the case. This is not surprising given the fact that the 
magnetic diffraction profiles are complex. Their shape is 
ascribed to the coexistence of long and short-range dy- 
namical correlations; see our discussion of Fig. [2] To the 
distribution of correlation lengths must correspond dis- 
tributions of spin-spin and spin-lattice relaxation rates. 
This may explain the exotic shape of the zero-field /xSR 
relaxation function. 

Relative to Tb2Sn207, the fluctuations probed by /xSR 



Paramagnetic correlations were studied by diffuse neu- 
tron scattering in Er2Ti2 07 crystals. Two scattering 
planes were investigated: (h, k, 0) and (h, k, k) at 2.00 (3) 
and 1.47(3) K, respectively. For the first (second) one 
a graphite (copper) monochromator was used delivering 
neutrons of wavelength 2.377 (1.275) A. No energy analy- 
sis of the scattered beam was performed. In order to deal 
with magnetic correlations only, additional maps were 
recorded at 50 K for both geometries and they were sub- 
tracted from the corresponding low temperature counter- 
parts. The resulting two maps were divided out by the 
square modulus of the Er 3+ magnetic form factor^ The 
effect of the form factor is modest: no more than 20% on 
the vast majority of the data. The resulting experimen- 
tal data are displayed in Fig. [14J It is to be noted that 
the intensities are negative for large regions in the two 
planes. This reflects the fact that the wavevector inde- 
pendent scattering associated with the CEF states of the 
Tb 3+ ions is larger at 50 K than at low temperature. 

Before attempting a quantitative analysis of the maps, 
a qualitative discussion is worthwhile. We first note the 
almost vanishing (2, 2, 2) spot. This means that the ferro- 
magnetic correlations are negligible. We observe a hexag- 
onal scattering loop in the plane (h, k, k) around the 
(2, 2, 2) position. Such a type of scattering is reminiscent 
of the intensity measured in the cubic spinel ZnC^CU 
(Ref.HH), in which the Cr ions also form a lattice of corner 
sharing tetrahedra. However, in the latter case the loop 
is in the plane (h, k, 0) and centered around (2, 2, 0). The 
scattering properties are therefore quite different for the 
two compounds. This reflects the difference in magnetic 
symmetry. The origin of the loops observed in ZnCr204, 
which were originally interpreted in terms of weakly in- 
teracting hexagonal spin clusters, is now taken as the sig- 
nature of extended exchange interactions for spin-ice and 
isotropic systems i 56 i In the following we show that the 
scattering loop in Er2Ti2 07 can basically be taken as a 
fingerprint of the properties of the exchange interactions 
within a single tetrahedron. 
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FIG. 14. (color online) Top two panels: magnetic diffuse neu- 
tron scattering intensity recorded for a crystal of Er2Ti207 
in the reciprocal (h, k, 0) and (h, k, k) planes at 2.00 (3) and 
1.47(3) K respectively. The positions in the reciprocal lat- 
tice are in 2n/a units, where a is the lattice parameter of the 
cubic unit cell. These maps are obtained as explained in the 
main text. To preserve the maps appearance, pixels with off 
scale intensities, e.g pixels influenced by Bragg reflections and 
critical scattering, as well as pixels located near the origin of 
the reciprocal lattice have been graphically eliminated: they 
are represented in white color. Bottom two panels: (h, k, 0) 
and (h, k, k) magnetic correlation maps computed with the 
tetrahedron model explained in the main text. The compari- 
son between the theoretical and experimental maps displayed 
above enables us to derive information on the Er2Ti2C>7 in- 
teraction constants. The lines drawn in the (h, k, 0) maps 
indicate the position of the cuts shown in Fig. 1161 

The discussion of the experimental results will be car- 
ried out in two steps. We shall first evaluate the mag- 
netic correlation length at the temperature of the mea- 
surements and then analyze the maps using a four-spin 
Hamiltonian. 
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FIG. 15. (color online). Magnetic scattering intensity ver- 
sus wavevector in the vicinity of the two reciprocal lattice 
positions qrjj^n where (h,k,l) = (1,1,1) and (2,2,0), respec- 
tively. The wavevector unit is 2n/a where a is the cube edge. 
The data are obtained from the intensity differences shown 
in Fig. Q3] by averaging the data at reciprocal space points 
located at equal distance from the two lattice positions. Two 
requirements were considered when choosing these two posi- 
tions: they do not lie close to the boundary of the maps and 
the magnetic intensity is relatively important. The lines are 
results from fits of Eq. [15] to the data. 

magnetic scattering curve is much larger than the instru- 
mental resolution, the convolution of Eq. [15] by the reso- 
lution function is unnecessary. The fits shown in Fig. 1151 
yield the magnetic correlation lengths £ m = k^ 1 = 3.6 (2) 
and 6.6(5) A for the (2,2,0) and (1,1,1) reflections mea- 
sured at 2.00 and 1.47 K, respectively. As expected, £ m 
shoots up as the sample is cooled toward the transition. 
These two values are comparable with the Er 3+ -Er 3+ ion 
distance d — 3.56 A. Hence the analysis of the experimen- 
tal maps shown in Fig. [TJ] can be performed considering 
the spin correlations within a single tetrahedron. This 
is the basis for our quantitative interpretation which is 
exposed below. 



1. Magnetic correlation length 

Here we determine the correlation length of the critical 
magnetic correlations. For this purpose we consider the 
scattered intensity measured in the vicinity of the recip- 
rocal positions q,(h,k,l) = 3(2,2,0) and 3(1,1,1) at T = 2.00 
and 1.47 K, respectively; see Fig. [TS] This critical scat- 
tering intensity is described by the sum of a Lorentzian 
function and a constant: 



£(|q- q(fc. 



kj) = 



1 



q(h,fe,ol 2 /«B 



h, (15) 



where K m is the inverse of the magnetic correlation 
length. The parameter Jl accounts for the magnitude 
of the Lorentzian, while Jo refers to a neutron inten- 
sity which is not related to critical scattering. Since at 
the temperature of experiments, the width of the critical 



2. Analysis of the diffuse scattering maps 

While our analysis of the magnetic scattering inten- 
sity in the vicinity of reciprocal lattice positions at low 
temperature shows that the measured wavevector depen- 
dence probes short-range correlations, a wavevector in- 
dependent scattering is also observed; see Fig. [15] This 
scattering reflects local physics, for example of crystal- 
field nature. Denoting M(q) a measured magnetic scat- 
tering map, we write 



M(q) =5 shlf t+S' Bca i e iV(q), 



(16) 



where S'shift accounts for the wavevector independent 
scattering, Sscaic gives the scale of the wavevector de- 
pendent magnetic intensity and N(q) is the prediction of 
our model that we describe now. 
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In Sec.|TI]we have mentioned that the Er 3+ crystal- held 
ground state doublet is well isolated since the first excited 
doublet is at about 74 K above the ground state in tem- 
perature units. Therefore for low temperature measure- 
ments, such as discussed here, it is a reasonable approx- 
imation to describe an Er 3+ ion as an effective S = 1/2 
spin. In Appendix lAl we discuss the physics of a tetrahe- 
dron of effective one-half spins embedded in a pyrochlore 
lattice. Assuming that only bilinear spin interactions are 
relevant, the Hamiltonian T-L t describing the interaction 
between the four effective spins can be written as a linear 
combination of four invariants; see Eq. lA38l Their weight 
is gauged by the four parameters Pi, i £ {1, 2, 3, 4}. 

The wavevector dependent scattering intensity result- 
ing from a single tetrahedron is proportional to^ 3 - 

AT(q)=5>^/^>/ m (q), (17) 
m 

where 

/m(q) = (18) 

a,f3 i,j n ^ "J 

The Debye- Waller factor is negligible for the region of 
the reciprocal space investigated and the temperatures 
at which the measurements were made. The symbols 
a, (5 refer to global Cartesian coordinates in the cubic 
unit cell. The relations between the local and global 
axes are explicitly given in Appendix IA II The indices 
i and j stand for the Er ions at the corners of the tetra- 
hedron, the positions of which are specified by Vi and 
Tj, and \n) and |m) refer to two of the sixteen single 
tetrahedron states, whose energy differences lie within 
the energy range across which the neutron scattering is 
integrated. The formula in Eq. [18] does not contain the 
Er 3+ form factor since it has already been divided out 
for the two presented maps. 

We have performed a global fit of the two recorded 
maps. It depends on two S s hm parameters, one per map, 
a unique scale parameter 5 sca i e and the four Hamiltonian 
parameters P{ . The best fit shown in Fig. [14] is achieved 
with the following interaction constants in kelvin units: 

Pi/fee = 10.4(6), P 2 /k B = -1.2 (5), 

P 3 /k B -3.0(3), P 4 /fc B = 8.4(1.3). (19) 

The quality of the fit can be assessed from cuts of the 
experimental and model maps as seen in Fig. 1161 

The fitting parameters described above were deter- 
mined from least square minimization, i.e. minimizing 
X 2 = [l/(N d - 7)] £(Mj - Ijf/rf with respect to the 
fitting parameters. Here aj is the statistical uncertainty 
on the neutron intensity Ij at the jth data point in the 
set of Nd data points and Mj is the model prediction (a 
function of the fitting parameters). An uncertainty S Pi 
on Pi is obtained by varying the neutron intensity by dlj 
at each of a large sample of the data points (one in eight) 
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FIG. 16. (color online) Cuts of the (h, k, 0) maps shown in 
Fig. [14] along the [1,1,0] and [0,1,0] directions. The circles 
represent the experimental data and the lines the results of 
the model. The origin of the horizontal axis is taken at the 
origin of the reciprocal lattice. 



and minimizing \ 2 m order to find the variation of each 
fitted parameter dPi. Then 8 Pi is estimated using the 

relation Sp = {l/N d ) ^^dP~JdI^. 

We have so far described the estimate of the statistical 
uncertainties. Three other origins for systematical uncer- 
tainties must be considered. First our model describes 
diffuse scattering and therefore pixels which are sizably 
influenced by critical scattering should be eliminated in 
the fitting procedure. For this purpose we have tested 
the influence of different radius cutoffs around the Bragg 
points, to the Pi values. This cutoff effect introduces the 
largest parameter uncertainties. Then the temperature 
at which a map is recorded is known with a finite un- 
certainty. Finally the uncertainty arising from the error 
bars on the matrix elements j'cef and icEF (see Eq. [2]) 
has also been assessed. The error bars given in Eq. [19] 
account for these statistical and systematical uncertain- 
ties. 

It is tempting to extract information from the Curie- 
Weiss constant using the high-temperature susceptibility 
formula derived Ross et aL— However, this should not 
be done since that formula depends on spectroscopic fac- 
tors which provide a description of the physics only at 
low temperature. A high-temperature expansion for the 
susceptibility is required to further analyze the data of 
Fig. [5] 

At this juncture it is of interest to mention a recent 
work by Savary et al. on Er2Ti20yi^ These authors an- 
alyze the spin wave dispersions measured for different 
orientations of the wave vector. The experiments were 
performed for temperature (30 mK) and field (3 T) val- 
ues where Er2Ti2C>7 is a polarized paramagnet. 20 The 
analysis of the data involves a nearest neighbor Hamil- 
tonian related to ours, with a different definition of the 
Hamiltonian parameters. Using the relations between the 
two sets of parameters given in Eq. IA411 and assuming 
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Pi = 2Pi we find for their parameters in kelvin units, 

Pi/fcB - 1-74 (1.26), P 2 /k B = -0.44 (74), 
P 3 /k B = 2.92 (34), P 4 /k B = 9.0 (1.0). (20) 

These P 2 , P3 and P4 values are in very good agreement 
with ours. However an important difference is observed 
in Pi, i.e. the parameter which controls the interaction 
along the local hard magnetic axis. We have simulated 
the diffuse scattering intensity obtained from this latter 
set of parameters and compared it with our experimental 
data. We obtain an acceptable fit of our data since the 
confidence parameter is % 2 = 1.44 to be compared with 
1.20 with the parameter set in Eq. Q15J At this point 
it must be noted that our data were recorded in zero 
external field, contrary to those of Ref. [26| and that the 
anisotropy ratio g±/g\\ = 4.3 that we have adopted is 
quite larger than the one chosen by Savary et ai: g±/g\\ 
= 2.4. 

As already mentioned the dispersion relation of the 
magnon modes has been computed in the framework of 
linear spin- wave theoryj2& Neglecting the gap, an assump- 
tion which is justified (see Sec. lIVB]) . we can compute the 
geometric mean v sw of the magnon velocities along the 
three directions of the Cartesian frame. The expression 
is: 

(Pi - Pi) 1 ' 2 (Pi - P 3 ) 1/6 (2Pi + P 3 ) 1/3 a 
Vsw 2 10 /3 3l/2 %■ 

(21) 

With the parameters given in Eq. [30] we find v sw — 
76(16) ms , in good agreement with the value v sw — 
84(2) and 82(2) ms -1 deduced in Sec. HVBl from the 
analysis of the specific heat data. 

Concerning the parameters derived from our diffuse 
scattering maps (Eq. [T9|) . the magnon velocity cannot 
be computed since (P4 — Pi) < 0. At first sight, it 
might suggest that the parameters we infer from the two 
neutron maps are not reliable. However, the expression 
written in Eq. [2U is deduced from a linear spin-wave 
approximation. This approximation might not be reli- 
able for a non-collinear magnet such as Er2Ti207. This 
argument is based on theoretical results for triangular 
magnets4^^60 



The critical exponent of the specific heat in the param- 
agnetic phase is typical for a three dimensional XY sys- 
tem. The specific heat data provides an estimate for the 
magnon velocity which is in accord with a recently pro- 
posed model based on linear spin-wave theory. Finally, 
concerning the bulk measurements, the magnetic phase 
diagram in the field-temperature plane has been deter- 
mined up to 1.7 T for the three main crystal directions 
of a cube. 

The /iSR data are consistent with the presence of a 
spin dynamics in the nanosecond time range in the or- 
dered state. It might be associated with the short-range 
correlations detected by neutron diffraction in addition 
to the long-range order =2 

Diffuse neutron scattering data recorded in the para- 
magnetic state were analyzed in terms of a Hamiltonian 
accounting for all bilinear interactions between the spins 
in a tetrahedron. The four symmetry allowed interac- 
tion constants were determined and three of them were 
found to be positive, with Pi and P4 being the largest. 
Note that the data analysis has assumed the interactions 
to be limited to the nearest neighbors Er 3+ ions. How- 
ever, it must be recalled that interactions between fur- 
ther neighbors might be important. For instance they are 
determinant for the type of magnetic order adopted by 
Gd2Ti2C>7 relative to Gd2Sn2 07 for which only nearest 
neighbors seem to matter. 9 The analysis of neutron scat- 
tering data for Dy2Ti2C>7 also confirms the importance 
of exchange interactions beyond nearest neighbors^ 

With the interaction constants and the spectroscopic 
factors known, it will be interesting to gauge any 
proposed theoretical phase diagram to the fact that 
Er2Ti2C>7 does order magnetically at Tn = 1.23(1) K. 
In addition, any reliable theory must be able to explain 
the short-range correlations observed below Tn and their 
nanosecond time scale dynamics. As it has been learned 
from the study of Tb 2 Sn 2 07, neutron-spin echo measure- 
ments might be useful to further characterize these exotic 
spin dynamics! 37 ' 53 The spin dynamics is the key feature 
which seems to distinguish a frustrated compound such 
as Er 2 Ti207 from a conventional magnet. 
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Appendix A: Physics of a tetrahedron of effective 
one-half spins embedded in a pyrochlore lattice 

Here we present a comprehensive quantum mechanical 
study of effective one-half spins embedded in a pyrochlore 
lattice. After specifying the geometry, we describe the in- 
variants which will enable us to build the Hamiltonian of 
the system. Then we explain the approximation that we 
use to compute the neutron diffuse scattering patterns. 
Finally, we compare our to others' equivalent Hamilto- 
nian operators. 



1. Geometry 

For completeness, we first provide a description of a 
tetrahedron in a pyrochlore lattice and of the spins at its 
corners. 

We choose a type A tetrahedron as the primitive unit 
cell; see Fig. [1] Each corner of the tetrahedron is occu- 
pied by a magnetic ion whose relative positions are given 
by n = f(0,0,0), r 2 = f(l,l,0), r 3 = f (1,0,1) and 
r4 = |(0,1,1), where a is the edge length of the cu- 
bic unit cell. The distance between two magnetic ions is 
d = a/(2\/2). At each corner of the tetrahedron there 
is a local D 3( i symmetry, where the C 3 axis is one of the 
cube diagonals. 

We denote the local C 3 axis at position i as Zi and 
define Xi and yi axes to form a local orthogonal basis. 
There is obviously some freedom in the choice of the Xi 
and yi axes. Below we list the unit vectors we have chosen 
for the bases at the four positions. For the first two 
positions we have 

xi = (l,l,-2)/V6, x 2 = (-l,-l,-2)/V6, 

yi = (-l,l,0)/V2, y 2 = (1,-1, 0)/A 

z 1 = (l,l,l)/V3, z 2 = (-l,-l,l)/V3, (Al) 

and at the remaining two positions we have 

x 3 = (-l,l,2)/V6, x 4 = (1,-1,2)/V6, 

y 3 = (l,l,0)/V2, y 4 = (-l,-l,0)/V2, 

z 3 = (-l,l,-l)/\/3, z 4 = (1,-1,-1)/V3. (A2) 

A given spin Si can be written using the cubic global 
axes as a basis or the local axes at position i. We use 
upper case superscripts to indicate components of the 
global basis, 

Si = Sf X + Sf Y + SfZ, (A3) 
and lower case subscripts for components of a local frame, 



Using the definitions above, we find 



Six — 


(Si + 


SY - 25f )/v / 6, 


(A5) 


Sly = 


1 nX 

{-st 


+ sY)/V2, 


(A6) 


Si z = 


/ n X 

(Si + 


S{ + S] f )/V3, 


(A7) 


&2 X — 


-(Si 


+ Sl +2S^)/V6, 


(A8) 


S 2y = 


(si - 


sl)/V2, 


(A9) 


S 2z = 


/ nX 

(-si 


- sl + s^)/V3, 


(A10) 


&3 X = 


(-sf 


+ s\ + 2si)/VE, 


(All) 


Szy = 


(S? + 


Sj)/V2, 


(A12) 


S 3z = 


( _oX 


+ sl - sD/Vs, 


(A13) 


Six = 




Sl + 2Si)/V6, 


(A14) 


Sly = 


( 


- Sj)/V2, 


(A15) 


Siz = 


(^4 


sj - sl)/Vz, 


(A16) 



and the inverse relations 



si 


= Si x /V6- 


Si y /V2 + Si z /V3, 


(A17) 


sl 


= Si x /Ve + 


Si y /V2 + Si z /V3, 


(A18) 


S Z i 


= -2S lx /V6 + Si z /V3, 


(A19) 


si 


= -S 2x /V& 


+ S 2y /V2-S 2z /V3, 


(A20) 


sl 


= -s 2x /Vg 


- S 2y /V2- S 2z /V3, 


(A21) 


sl 


= -2S 2x /V6 + S 2z /V3, 


(A22) 


si 


= -s 3x /Ve 


+ S 3y /V2-S 3z /V3, 


(A23) 


oY 
s 3 


= s 3x /Ve + 


S 3y /V2 + S 3z /V3, 


(A24) 


oZ 
°3 


= 2S 3x /V6- 


- S 3z /V3, 


(A25) 


aX 
04 


= Six/vQ- 


S4y/V2 + S/1 Z /V3, 


(A26) 


S4 


= -s^x/Vq 


- Siy j V2 — Si Z j VO. 


(A27) 


S4 


= 2S ix /V6- 


- S iz /V3 


(A28) 



2. Determination of the invariants 



The interaction Hamiltonian T-L between the spins in 
the lattice is described by bilinear operators of the form 
SiaSjp where i and j are nearest neighbors. The space 
group Fd3m allows four different terms for a Kramers 

ion 61-63 



Si — Si x ii + SiyYt + Si Z Zi. (A4) 



H = P1X1 + P 2 X 2 + P 3 X 3 + P 4 X A , (A29) 
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where Pi are constants and 



3. Single tetrahedron approximation 



X 1 = ~J2S iz S j „ (A30) 

<%2 = g- [^4j(SizSj+ + Sj z Si+) 

+A* j (S iz S j - + S jt Si-)] , (A31) 

X ^ = \J2 ( 5 *+ S J+ + A v Si-Sj-), (A32) 



(A33) 



The sum is over all the pairs of nearest neighbors (see 
Eq. [S] for an example of the use of the notation in the 
case of the Heisenberg interaction). Below we give the 
expressions of the phase factors Ay introduced for X% 
and A3: 



A12 = A 34 = 1, 

A13 = A 24 = e = exp 

A14 = A 23 = e* = exp 



2ni 

TJ' 

47TI 

~3~ 



(A34) 
(A35) 

(A36) 



Note that the sum of all four invariants is the isotropic 
exchange interaction YU=i X = Yllij) ' &j (obtained 
when P 1 = P 2 =P 3 = P 4 ). 

To get some insight in the % expression, as an example, 
we consider the product Si • S2. Let us first focus on the 
contribution of the first invariant to this product: 

X\ : ^"g - ^ 2 (^ lz ^ 2z ^2zSix) = ——Sx z S2z- 
The contribution of the second invariant is easily found: 

A2 : —2-^- (Si z S2x + S2zS\ x ) . 

In the same way the contributions of A 3 and A4 can be 
derived. They are written in terms of the local Carte- 
sian axes. Transforming to the global axes we derive the 
expected relation 

So qX c?X 1 nY qY 1 qZ qZ 

1 ' ^2 — *^2 ' ^1 ^2 ' ^1 ^2 • 



In its most general form (Eq.[A29| H includes symmet- 
ric (Heisenberg) and the antisymmetric (Dzyaloshinskii- 
Moriya) interactions. For completeness, we give the ex- 
pression of the Dzyaloshinskii-Moriya Hamiltonian. In 
terms of the invariants, we derive 



"Hdm — — -St, dm ( 4Afi — -A 2 + A 3 — 2X4 J , (A37) 

where i?t,DM scales the Dzyaloshinskii-Moriya interac- 
tion. 



In order to compute diffuse neutron scattering pat- 
terns exact eigenstates of H must be used, however, 
the Hamiltonian given by Eq. IA29I is unsolvable in gen- 
eral. Therefore instead of using the full Hamiltonian 
(Eq. IA29[) . we find exact solutions to an Hamiltonian 
Ht restricted to tetrahedra of a single type. This ap- 
proach, which has been used to analyze diffuse neutron 
scattering patterns ) 61,64 is a considerable simplification 
of the original model. It amounts to replacing the sums 
in Eq. IA29I over all tetrahedra, as it is implicit from the 
definition of the x-j's, to a sum over all A- type (or all B- 
type) tetrahedra (Fig. [T]), as explained in Ref. |62j. Thus 
each only three of the six nearest neighbors of each Er 
atom are included in the calculation. Instead of the orig- 
inal interaction constants Pi we introduce a similar set 
of constants Pi, such that 



n t = 



E 

A tetrahedra 



P1A1 + P2A2+P3A3 + P4A4. (A38) 



Here the sums appearing in the Xi operators are limited 
to nearest neighbor spins belonging to single A tetrahe- 
dra, unlike in Eqs. lA30llA33l While the full ramifications 
of this approximation are not understood, at the very 
least we can estimate that the interaction constants P; 
are approximately a factor of two larger than Pi to com- 
pensate for the missing exchange paths in Eq. IA38I 



4. Relations between Hamiltonian parameters 

Other groups have proposed an Hamiltonian for the 
description of an effective one-half spin system such as 
Er2Ti2C>7 for which the rare-earth crystal-field ground 
state is a Kramer's doublet. 

Ross et al£& and Savary et alM^ write the Hamilto- 
nian as 



u 



{J zzSizS j z J± {Si+Sj— + Si—Sj+) 



+ J±± (lijSi+Sj+ + l*jSi-Sj-) 

+ J z ± [Su + QSjJ) (A39) 

where 7 and £ are 4x4 matrices: 



( 


-1 


e i7r/3 g-i-7r/3\ 


-1 





g-i7r/3 gi^/3 


e i7r/3 


e — t7r/3 


-1 


\ e -w 3 


e w/3 


-1 J 



(A40) 

Taking into account (i) the different labeling in the 
four sites of the tetrahedron as well as (ii) the different 
choice for the definition of the axes perpendicular to the 
local threefold axis which are adopted in these references 
compared to ours, the relations between the Hamiltonian 
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parameters are 



■hz — — g-Pl) Jz± — —P2, 

J±± = \h, J±= \h- (A41) 

As explained in Sec. IA 31 we expect Pi s» 2iV 

Considering Yb2Ti207 and others crystal-field ground- 
state doublet pyrochlore compounds, Onoda and 
Tanaka££~— have also used an Hamiltonian equivalent to 
Eq. IA391 but with slight differences in the Hamiltonian 
parameters labeling. 

Two other research groups have studied the interac- 
tion of spins on the pyrochlore lattice with the purpose 
to extract values of the interaction constants. However, 
since their interest was on the analysis of relatively high- 
temperature data, the full angular-momentum was 
used] 69 ! 70 Here we have recorded the diffuse scattering 
intensity at a sufficiently small temperature that it is 
justified to work within the effective one-half spin frame- 
work. It is possible to project the full angular momen- 
tum into the ground-state Kramer's doublet. However, it 
seems that only semiformal relations between the param- 
eters of the effective one-half spin and the full angular- 
momentum models can be derived^ Therefore we do not 
consider these relations in our analysis in Sec. IV B 21 



Appendix B: Estimate of the electric field gradient 
acting on the 1B7 Er nucleus 

Here we provide an estimate for the principal value of 
electric field gradient tensor V zz . It is required for the 
computation of the nuclear contribution to the specific 
heat. 

In an insulator V zz is written as the sum of two terms 
V^ tt and V Z J , respectively, modeling the lattice and 
4/-shell contributions. The former contribution is ex- 
pressed as V zz u = — ■ \~ ^ where AP> is a crystal- 
electric field parameter and 7oo and ui are a Sternheimer 
coefficient and the screening coefficient of the crystal 
field, respectively. From the literature values A = 
41.5 (1.1) meVao 2 (Ref.[H) where a = 52.92 pm is the 
Bohr radius and (l-^ 00 )/(l-a 2 ) = 210 (30) (Ref.©, we 
obtain V^f* = -1.24 ( 21) x 10 22 Vm" 2 . The 4/-shell con- 
tribution is written V Z J = —4^^2(1 — Ro){f'~ 3 )ifE i f. 
Rq is a Sternheimer coefficient, 62 is a Stevens coeffi- 
cient, and (r~ 3 )4f and E±f are the expectation values 
respectively of the cube of the inverse distance between 
the nucleus and the Er 3+ 4/ shell, and the quadrupole 
operator 3 J 2 — J (J + 1) acting on the 4/ shell. As usual, 
£0 is the permittivity of free space. From the literature 
values R Q = 0.29(1) (Ref.[nl), d 2 = 2.54 x 10" 3 , (r' 3 )^ 
= 11.36 3 (Ref.[zi) and E 4f = -7.86 (Ref.[H), we get 
V*/ = 1.6 x 10 21 Vm~ 2 . Summing up the two contribu- 
tions we obtain V zz = -1.1 (3) x 10 22 Vm~ 2 . 
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